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o : 

We study the quark-gluon mixed condensate gl^qa^G^q), using the SU(3) C lattice QCD 
$H ' with the Kogut-Susskind fermion at the quenched level. We generate 100 gauge configura- 

tions on the 16 4 lattice with j3 = 6.0, and perform the measurement of giqa^G^q) at 16 
points in each gauge configuration for each current quark mass of m q = 21,36,52 MeV. 
Using the 1600 data for each m q , we find rriQ = giqu^G^q) / (qq) ~ 2.5 GeV 2 at the lattice 
scale of a -1 ~ 2GeV in the chiral limit. The large value of g^qa^G^q) suggests its impor- 
tance in the operator product expansion in QCD. We study also chiral restoration at finite 
temperature in terms of g^qa^G^q) , which is another chiral order parameter. We present 
t-H ' the lattice QCD results of giqcr^G^q) at finite temperature. 

> : 

§1. The importance of the quark-gluon mixed condensate g (qcr ^ U G 

cn \ 

^) \ The non-perturbative nature of quantum chromodynamics (QCD) is character- 

ized by its nontrivial vacuum structure such as various condensates. Among various 
condensates, we emphasize here the importance of the quark-gluon mixed conden- 
sate g(qo ^ U G ' ^q) = g(qo~^uG^ u ^X A q). First, the mixed condensate represents a 

,— j H ■ direct correlation between quarks and gluons in the QCD vacuum. In this point, 

the mixed condensate differs from (qq) and (G^G^) even at the qualitative level. 
Second, this mixed condensate is another chiral order parameter of the second lowest 
dimension and it flips the chirality of the quark as 

• i— i . 

giqo-^G^q) = g(q R (a^G^) q L ) + g(q L (a^G^) q R ). (1-1) 

Third, the mixed condensate plays an important role in various QCD sum rules, 
especially in the baryonsP'^ the light-heavy mesonP^ and the exotic mesons.^ In 
the QCD sum rules, the value m\ = g(qcr ^ V G ^ v q) / (qq) ~ 0.8 ± 0.2 GeV 2 has been 
proposed as a result of the phenomenological analyses EJ^EMH) However, in spite of 
the importance of g(qc r ^G flu q}, there was only one preliminary lattice QCD studyP 
which was performed with insufficient statistics (only 5 data) using a small (8 4 ) and 
coarse lattice (/3 = 5.7). 

Therefore, we present the calculation of g(qcr^ v G^ u q) in lattice QCD with a 
larger (16 4 ) and finer (j3 = 6.0) lattice and with high statistics (1600 data). We 
perform the measurement of g(qa ^ U G ^q) as well as (qq) in the SU(3) C lattice at the 
quenched level, using the Kogut-Susskind (KS) fermion to respect chiral symmetry. 
We generate 100 gauge configurations and pick up 16 space-time points for each 
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configuration to calculate the condensates. With this high statistics of 1600 data for 
each quark mass, we perform a reliable estimate for the ratio ra\ = g(qo ^G ^q) / (qq) 
at the lattice scale of a -1 ~ 2GeV in the chiral limit 

§2. Lattice Formalism 

Since both of the condensates (qq) and g(qo ^q) are chiral order parameters, 
they are sensitive to explicit chiral-symmetry breaking. Therefore, we adopt the KS 
fermion, which preserves the explicit chiral symmetry in massless quark limit, m = 0. 
On the other hand, the Wilson and the clover fermions would not appropriate for 
our study, because these fermions explicitly break chiral symmetry even for m = 0. 

The action for the KS fermion is described by spinless Grassmann fields X , x an d 
the gauge link- variable = exp[— iagA^]. When the gauge field is set to be zero, 
the SU(4) f quark-spinor field q with spinor i and flavor / is expressed by x as 



ll(x) 



lY.( r phx(x + P ), r p = 7 f 7 f 7 f 3 7f, p=(pi,P2, P 3,PA) (2-1) 



where p with p^ G {0, 1} runs over the 16 sites in the 2 4 hypercube. When the gluon 
field is turned on, we insert additional link-variables in Eq. (|2-lj) in order to respect 
the gauge covariance. Hence, the flavor-averaged condensates are expressed as 

a 3 (^) = -^Tr[(^(x)^(x)}] =-l^Tr[r p r p t ( x (x + p) X (x + p))] ,(2-2) 
/ P 

a 5 g(qcrf,uG^ u q) = -- EE Tr ( x )) a ^ G ^ 

f W 

= E E Tr [ U ±nM* + P) r p' r l U(x + p'Mx + p)) <V Gf v (x + p) ] , 
pv p 

(2-3) 

p' = p±p±u, 

where the sign ± is taken such that the sink point (x + p') = (x + p± p±v) belongs 
to the same hypercube of the source point (x + p). Here, in order to respect the 
gauge covariance, we have used, in Eq. ()2-3|) . 

U±^ ±u (x) = ~ [ U±^x)U ±u (x ± M ) + U ±v (x)U ±fl (x ± v) ] , (2-4) 

where the definition of U-^(x) = Up[(x — p) is used. 

On the gluon field strength G^, we adopt the clover- type definition on the 
lattice, 

G>(«) = ^ E XA Tr i xA i U ^( s ) + + U-p-u(s) + U^,(s)} - X A {p <-> u}] , 
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which contains no 0(a) discretization error. This is also an advanced point which is 
absent in Ref . EJ) . 

§3. The lattice QCD results 



We calculate the condensates (qq) and g(qa ^G ^q) using the SU(3) C lattice 
QCD at the quenched level. We perform the Monte Carlo simulation with the 
standard Wilson action at f3 = 6.0 on the 16 4 lattice. The lattice unit a ~ 0.10 fm is 
obtained so as to reproduce the string tension a = 0.89 GeV/fmPJ We use the quark 
mass m = 21,36,52 MeV, i.e., ma = 0.0105, 0.0184, 0.0263. For the fields x, X, 
the anti-periodic condition is imposed. We measure the condensates on 16 different 
physical space-time points x in each configuration as x = (x\, X2, x%, X4) with x^ € 
{0,8} in the lattice unit. For each m, we calculate the flavor-averaged condensates, 
and average them over the 16 space-time points and 100 gauge configurations. 

Figure n shows the values of the bare condensates a 3 (qq) and a^g(qa^ v G^ v q) 
against the quark mass ma. We emphasize that the jackknife errors are almost 
negligible, due to the high statistics of 1600 data for each quark mass. From Fig. ^ 
both (qq) and g(qa ^G '^q) show a clear linear behavior against the quark mass m. 
Therefore, we fit the data with a linear function and determine the condensates in 
the chiral limit. The results are summarized in Table [IJ 
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Fig. 1. The bare condensates a 3 (qq) and a J g(qa IJ , v G ^q) plotted against the quark mass ma. The 
dashed lines denote the best linear extrapolations, and the cross symbols correspond to the 
values in the chiral limit. The jackknife errors are hidden in the dots. 



Table I. The numerical results of a (qq) and c 1 g(qo ^q) for various ma. The last column 
denotes their values in the chiral limit obtained by the linear chiral extrapolation. 



ma = 0.0263 


ma = 0.0184 


ma = 0.0105 


chiral limit 


a^qq) -0.04240(16) 
a 5 g(qay,„Gy.vq) -0.01882(15) 


-0.03247(15) 
-0.01498(14) 


-0.02212(16) 
-0.01088(14) 


-0.00872(17) 
-0.00565(14) 



We check the reliability of our lattice QCD results by considering the finite 
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volume artifact. In order to estimate this artifact, we carry out the same calculation 
imposing the periodic boundary condition on the Grassmann fields % and x, instead 
of the anti-periodic boundary condition, keeping the other parameters same. We 
obtain that the results with different boundary conditions almost coincide within 
about 1% difference, and thus we conclude that the physical volume V ~ (1.6 fm) 4 
in this simulation is large enough to avoid the finite volume artifact.^ 

The values of the condensates in the continuum limit are to be obtained after the 
renormalization, which, however, suffers from the uncertainty of the non-perturbative 
effect. As a more reliable quantity, we provide the ratio m\ = g{qo' j j ai G lxv q) / '{qq), 
which is free from the uncertainty from the wave function renormalization of the 
quark. 

Now, we present the estimate of m\ using the bare results in SU(3) C lattice QCD 

as 

ml = g{qa^G^q)/{qq) ~ 2.5 GeV 2 (/? = 6.0 or a' 1 ~ 2GeV). (3-1) 

The large value of m\ suggests the importance of the mixed condensate in OPE. 
Although we do not include renormalization effect, this result itself is determined 
very precisely.^ 

§4. Discussions and Outlook 

For comparison with the standard value in the QCD sum rule, we change the 
renormalization point from /i ~ 7r/a to fi ~ 1 GeV corresponding to the QCD sum 
rule. Following Ref. EJ, we first take the lattice results of the condensates as the 
starting point of the flow, then rescale the condensates perturbatively. We adopt the 
anomalous dimensions at the one-loop level and choose the parameters ^Iqcd = 
200 — 300MeV and Nf = corresponding to quenched lattice QCD. We obtain 
m o|^=iGeV = 9{q^^uG^q)/{qq)\^ =lGeW ~ 3.5 - 3.7 GeV 2 . Comparing with the 
standard value of m\ = 0.8 ± 0.2 GeV 2 in the QCD sum rule, our calculation results 
in a rather large value. (Note that the instanton model has made a slightly larger 
estimate as m 2 , ~ 1.4 GeV 2 at fj, ~ 0.6 GeV.^) For a more definite determination 
of ?71q, the non-perturbative renormalization scheme may be desired. 

We again emphasize that the mixed condensate giqcr^G^q) plays the important 
roles in various contexts in quark hadron physics. Hence, it is preferable to perform 
further studies. In particular, the thermal effects are interesting in relation to chiral 
restoration, because the mixed condensate is another chiral order parameter. Con- 
sidering also the on-going experiments for the finite-temperature QCD in the RHIC 
project, we investigate the thermal effects on the mixed condensate g(qcr^ u G^ u q). 
We perform the calculation using the 16 3 x N t lattices with N t = 16, 12, 10, 8, 6, 4 at 
P = 6.0. Figure |21 shows the lattice QCD results for the mixed condensate at finite 
temperature. We find a drastic change of the mixed condensate around the critical 
temperature T c , which reflects chiral-symmetry restoration.^ 

In summary, we have studied the quark-gluon mixed condensate giqcr^G^q) 
using SU(3) C lattice QCD with the KS fermion at the quenched level. For each quark 
mass of m q = 21,36,52 MeV, we have generated 100 gauge configurations on the 
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Fig. 2. The quark-gluon mixed condensate g{qo ^vG llv q)T /g{qa^G,j,uq}T^o plotted against the 
temperature T. The jackknife errors are hidden in the triangles. The vertical dashed line 
denotes the critical temperature T c ~ 280MeV at the quenched level. 



16 4 lattice with = 6.0. Using the 1600 data for each m g , we have found mg = 
g{qo~^uG^ v q) I (qq) ~ 2.5 GeV 2 in the chiral limit at the lattice scale corresponding to 
(3 = 6.0 or a -1 ~ 2 GeV. We have also investigated g{qa^ u Gfj, u q) at finite temperature 
in lattice QCD. 
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